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ABSTRACT 



Oscillations were observed across the whole solar disk using the Doppler shift 
and line depth of spectral lines from the CO molecule near 4666 nm with the 
National Solar Observatory's McMath/Pierce solar telescope. Power, coherence, 
and phase spectra were examined, and diagnostic diagrams reveal power ridges 
at the solar global mode frequencies to show that these oscillations are solar 
p-modes. The phase was used to determine the height of formation of the CO 
lin es by comparison with the IR continuum intensity phase shifts as measured 



m 



Kopp et al.l (119921 ): we find the CO line formation height varies from 425 < 
z < 560 km as we move from disk center towards the solar limb 1.0 > fi > 0.5. 
The velocity power spectra show that while the sum of the background and p- 
mode power increases with height in the solar atmosphere as seen in previous 
work, the power in the p-modes only (background subtracted) decreases with 
height, consistent with evanescent waves. The CO line depth weakens in regions 
of stronger magnetic fields, as does the p-mode oscillation power. Across most of 
the solar surface the phase shift is larger than the expected value of 90 degrees 
for an adiabatic atmosphere. We fit the phase spectra at different disk positions 
with a simple atmospheric model to determine that the acoustic cutoff frequency 
is about 4.5 mHz with only small variations, but that the thermal relaxation 
frequency drops significantly from 2.7 to mHz at these heights in the solar 
atmosphere. 
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Introduction 



The solar infrared absorption lines from the fundamental band of the CO molecule at 
4666 nm are well-known diagnostics of the temperature minimum and chromosphere of the 
solar atmosphere. I nitial discovery of off-limb emission in the quiet Sun from these lines 
( iNoyes &: Halllll972[) instigated the deve l opment of thermally bifurcated models of the so- 



lar atmosphere, ( IWiedemann et al.lll994j : lAyres fc Rabinlll996l ) and the latest models draw 



(Ascensio Ramos et al. 2003: Wedemever-Bohm et al. 


2005 


). The CO lines are known to 


show 5 minute period Doppler oscillations ( 


Noves & Hall 


1972: Uitenbroek, Noves & Rabin 


1994: Solanki et al. 1996: Uitenbroek 2000c 


i) which are presumably driven from the solar 



global p-modes, and images of the solar surface in the lines r eveal that rnagnet ic regions of 



the solar atmospher e show a reduced m olecular line strength (jUitenbroekll2000al ). Models of 



the line formation ( lUitenbroekl l2000al ) suggest the formation height of the strong CO lines 
at 4666 nm, particularly the 3-2 R14 transition at 4665.8 nm, varies from roughly 150 to 
580 km above the t^qq = 1 formation height (z=0). 

Studies of the solar p-modes at various heights in the solar atmosphere have been done 
using a variety of spectral lines. It is thought that the p-modes are reflected back into the 
Sun by the steep gradient of density at the photosphere, and the wave solution requires 
an evanescent component (characterized by an exponential amplitude drop and a constant 
velocity phase with height) to be present above the reflecting layer. Simultaneous mea- 
surements in various spectral lines have shown that for the solar p-mode fr equencies the 
velocity oscillations are in phase at different heights (iLites fc Chipmanlll979[) but that the 
amplitude of the velocity oscillation increases w ith height ( ICram. Brown &: Beckers! 119771 : 
Ruiz Cobo. Rodriguez Hidalgo fc Colladoslll997| ) 



In addition to velocity oscillati ons, intensi t y osc illations are also observed in the solar 
atmosphere; since the early work of ISchmeiden (119761 ) the relationship between the two has 
become an important diagnostic of the physical conditions of the solar atmosphere. The 
phase between the temperature and velocity oscillations at various heights in the solar at- 
mosphere has been probed by many authors thro ugh the use of line depth and Doppler shift 
measurements of a different solar spectral line s (jPeubner et al.l Il992l : lAl. Bendlin fc Kneer 



19981 : lOliviero et al.l Il999l : IStrauss et al.lll999[ ). A motivation here is to understand if the 



solar atmosphere reacts adiabatically to the p-mode oscillations or not, and if not, then to 
determine what other physical processes may be important. Since the velocity oscillation am- 
plitudes increase with height, the solar atmosphere undergoes compression when the velocity 
oscillations move the gas downwards towards the photosphere, as the higher layers catch up 
to lower layers. If energy is conserved by the atmosphere then the gas temperature will reach 
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a maximum value when the compression reaches a maximum; this corresponds to the time 
when the gas velocity is zero just before it begins to move away from the solar photosphere. 
If the gas temperature is plotted versus the atmospheric velocity (defined as positive when 
the gas moves away from the photosphere, which is the opposite direction as the observed 
Doppler shift) in an adiabatic atmosphere there will be a 90 degree phase shift between the 
two quantities, with the temperature peak leading the atmospheric velocity peak. The phase 
between the temperature and the atmospheric velocity is referred to as the "I-V phase". In 
summary, observations from many different spectral lines show that the I-V phase is less 
than 90 degrees in the photosphere, around 90 degrees from the temperatu re minimum to 



the middle chromosphere, and above 90 de grees in the high chromosphere (jPeubner et al. 



1990l : iMassiello. Marmohno fc Strausslll998h 



Complicating the spectral line measurements of I-V phase is the fact that chromospheric 
spectral lines are often formed in a non-LTE environment, and the spectral line depth can 
change non-linearly as the tempera t ure of the solar at mosphere changes. Investigations from 



Kopp and co-workers (IKoppI Il990l : iKopp et al.l [1992) addressed this problem by observing 



temperature oscillations using thermal infrared continuum channels at wavelengths from 
50 to 800 microns from the airborne NASA Kuiper Astronomical Observatory (KAO). The 
continuum intensity depended only on the continuum opacity, and so the observed intensities 
were linearly related to the temperature of the solar atmosphere and provided excellent 
thermal probes. These wavelengths sampled the solar atmosphere from heights of roughly 340 
to 820 km above z=0 and thus the temperature minimum through low chromosphere could 
be explored. Since the continuum observations lacked a way to determine the gas velocity, 
simultaneous ground-based spectral observations were made using both photospheric and 
chromospheric spectral lines. 

The KAO observations showed an I-V phase shift consistent with results from spectral 
lines which sampled the same atmospheric heights. In particular, the I-V phase was observed 
to be 122 degrees at 340 km, and decreased with height to approach 95 degrees at 600 km; 
the solar atmosphere became more adiabatic as one moved from the temperatu re minimum 



up to the lower chromosphere. In order to explain the non-adiabatic quality, iKopp et al. 



( 119921 ) introduced a simple thermal relaxation term into the hydrostatic gas thermal balance 
equation and solved for the phase in a simple model assuming an isothermal plane-parallel 
non-magnetic atmosphere. The resulting analytic expressions were used to determine the 
values of the thermal relaxation frequency ujr- 

In the following paper we discuss our observations of temperature and Doppler oscilla- 
tions using the 4666 nm absorption lines of CO. Since the CO spectral lines are formed in 
LTE conditions, the line depth provides a more direct probe of the gas temperature than 
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other chromospheric lines; and as spectral lines their Doppler shift gives a direct measure- 
ment of the atmospheric velocity and thus they are superior to observations using infrared 
continuum channels. In Section 2 we discuss the data collection, Section 3 we outline the 
data reduction and analysis procedures, and in Section 4 we discuss the results obtained from 
the intensity, velocity, and I-V phase measurements from these lines as they vary across the 
solar surface and with height in the solar atmosphere. We apply the simple analytic model 
of Kopp to our I-V phase spectra to estimate the gas properties in the solar atmosphere. 



2. Observations 

The observations presented here were collect ed on 16 Sep tember 2009 at the National 



Solar Observatory McMath-Pierce solar telescope (jPiercelll964l ) using the 13.5 m main spec- 
trograph. In order to image the full solar disk onto the detector, the "integrated light" flat 
mirror was inserted into the normal optical path of the telescope, thus providing unfocused 
sunlight into the observing room. The solar image was formed with a 70 mm diameter 
1400 mm focal length (effective focal length at 4666 nm) CaF2 singlet lens, which formed 
a 13.2 mm diameter prime focus image on the slit of the main spectrograph. The resulting 
f/20 beam overfills the spectrograph coUimating mirror, reducing the transmitted solar flux 
and increasing the spectrograph stray light, since the spectrograph is designed to work at 
f/60. These sacrifices were made in order to quickly scan the entire solar disk. 

At the spectrograph exit port, the NSO Array Camera (NAG) was fed with a single fold 
mirror. The NAG is a cryogenically cooled Aladdin III InSb 1024x1024 pixel array detector. 
The beam is fed through a cold (55K) order sorting filter to a single cold fold mirror and then 
to the InSb array; for these observations a 25 mm diameter filter centered at 4666 nm with a 
band-pass of 17 nm was used. The NAG was operated in "streaming" mode where about 4.5 
frames per second were collected and stored, during which time the solar image was scanned 
back and forth across the spectrograph slit asynchronously by moving the McM-P heliostat 
mirror (Ml). The scanning sequence took about 50 seconds to cross the solar disk, and 160 
scans of the solar disk were taken during 133 minutes from 1647 to 1900 UT. 

The GaF2 lens determines the spatial sampling along the spectrograph slit, and here we 
obtain 3.92 arcsec per pixel. The scan rate and camera readout time determine the spatial 
sampling in the direction perpendicular to the spectrograph slit, and there we achieve 12.39 
arcsec per pixel. At this wavelength the spatial resolution of the 70 mm diameter lens is 
13.75 arcsec and so the image is oversampled along the slit but not Nyquist sampled in the 
scan direction. The dispersion of the spectrograph is 0.0067 nm per pixel, but the spectral 
resolution is estimated to be ten times less at i? = 60, 000 by comparing the observed spectra 
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with model spectra provided by Uitenbroek. Simulations degrading the model spectra to this 
resolution also show that there seems to be significant scattered light, likely caused by feeding 
the spectrograph with a faster beam than its design and the resulting internal reflections. 

3. Data Reduction and Analysis 

The NAC was used in single-frame read streaming mode instead of using double- 
correlated sampling; in this way the InSb array pixels were read-out just once per exposure 
and the cadence was fast. NAC exposure times were very short, only 50 milliseconds per 
frame. Tests of this mode have determined that the entire array bias level changes slightly 
from frame to frame, but that individual pixel biases remain constant. Thus in addition to 
the usual dark current subtraction procedure, a set of pixels which were not exposed to light 
from the spectrograph were used to remove the changing bias level. A list of bad pixels (less 
than 300 of the one million pixels) was computed from a dark frame and the values of these 
pixels were replaced with a linear extrapolation from neighboring pixels. Gain correction 
was simply accomplished by scanning the spectrograph grating while taking a sequence of 
images, and selecting the maximum intensity value for each pixel out of the resulting data 
cube. This intensity level occurs when a pixel is exposed to the local continuum, and in 
this way the spectral lines are removed from the gain correction image. The slowly changing 
transmission profile of the cold order-sorting filter remains in this procedure, but it is re- 
moved with a linear fit to the intensity of the pixels which observe the continuum spectrum. 
Figure 1 shows a sample spectral frame with several spectral lines labeled. The signal-to- 
noise ratio of the data can be estimated with the standard deviation of the continuum pixel 
intensities to be roughly 50. 

Three strong CO spectral lines were then fit in the corrected data frames. Routines 
were used to fit both Gaussian (from IDL) and Voigt functions (written in IDL by D.M. 
Zarro) to the fine profiles. Initial guess for the position of each spectral line accounted for 
the spectrograph drift seen in the data (roughly 0.01 nm per hour) and the Doppler shift 
from solar rotation; in this way neighboring solar or atmospheric spectral lines were avoided. 
As part of the Voigt fitting, the parameters from the Gaussian fit were used as a first guess. 
Because of this the Voigt fitting procedure was more sensitive to noise. A comparison of 
the line center wavelengths returned from the Gaussian and the Voigt fitting was done; no 
systematic differences were seen between the values, and the small differences between the 
line centers were consistent with expectations from spectral measurements with a signal 
to noise of about 50. Since the simpler Gaussian fitting was able to fit more spectral lines 
successfully and no differences were seen between the two methods when both fits converged. 
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the Gaussian fit values were used in the rest of the analysis. 



Scanning with the heliostat mirror (Ml) at the McM-P without any limb guiding in- 
troduced significant image drift (about 100 arcsec E/W and 20 arcsec N/S) during the 2 
hour observing run. Because the image scanning extended above the solar limbs, the image 
drift was accommodated within the field-of-view of the observation. The solar images were 
registered using the continuum solar disk. Because the massive 2 m diameter Ml has a 
limited scan speed, the image scanning proceeded East to West in one scan, and then back 
West to East in the following scan, rather than having a "flyback" to scan continuously in 
one direction. In this way pixels on the solar disk are sampled at different times. To correct 
for this, alternate pairs of scanned images were linearly interpolated to form a new image 
at a single mean time. Finally, irregularities and dust on the spectrograph slit produced a 
streaking in the data, and this was removed. A co-registered and calibrated time series of 
the spectral hne depth, the hne position and line width was produced for each spectral line. 
In order to improve the signal to noise ratio of these parameters, the quantities from all 
three strong CO lines were averaged together. While the lines may sample slightly different 
regions of the so lar atmosphere, since the d epths of the three lines vary by less than 3% 
from each other (jWallace fc LivingstonI 120031 1 . we assume that they originate from roughly 
the same heights. 



Following iKopp et al.l ( 119921 ) a set of 10 pixels at disk center were averaged together, 
and the mean atmospheric velocity and line depth were plotted as a function of time. A 
cross-correlation function was computed between the velocity and line depth; the peak of the 
cross-correlation function revealed that the temperature led the atmospheric velocity signal 
on average by a little more than 90 degrees (assuming a dominant oscillation period of 5 
minutes). This cross-correlation procedure was then repeated for each pixel across the solar 
surface, and a map of the resulting I-V phase was produced. 

A more detailed analysis of the I-V phase was produced with a Fourier decomposi- 
tion of the velocity and line dep t h tirn e ser ies this analy s is fol lows standard techniques 
used by iLites. Rutten &: KalkofenI (119931 ) and iKrijger et al.l (120011 ) . The individual time se- 
ries were transformed to produce maps of the power in the atmospheric velocity and line 
depth parameters across the solar disk. If /^(t) is the time series of the atmospheric veloc- 
ity, we define F^,(i') as the Fourier transform of the atmospheric velocity. As discussed in 
Krijger et al.l (120011 ). we can decompose the power spectrum into real and imaginary com- 
ponents F^i^u) = Qyi^i') + ihy{i'). We define similar quantities for the line depth as fdit) and 
Fdiv). The cross-power spectrum was computed using CPdviy) = Fd{i^)F*{h'), where F*^^) 
is the complex conjugate of the transform of the velocity time series. 



From these we now compute the coherence and the phase spectra in the normal way, ex- 
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cept that the particular way of computing the ensemble averages deserves mention. We can 
express the cross-power spectrum in terms of real and imaginary components as CPdv{i^) = 
Cdv{^) + iddv{i^)- We use these complex coefficients for the ensemble averaging. In par- 
ticular, for the center-to-limb phase analysis, we average the Cdv{i') and ddv{^) coefficient 
in each bin of pixels at a given value of /i, and then compute the phase spectrum using 

= atan(^^^^). For the diagnostic diagrams we average the complex coefficients over 
and ky (the wavenumbers in the N/S and E/W solar directions) durin g the comput a tion o f 
the average horizontal wavenumber kh = ^/^^ + k^. As pointed out in lKrijger et al.l ( I2OOII ). 
averaging the complex coefficients before computing the phase reduces the noise in the cal- 
culation, whereas computing the phase first and averaging the noisy phase values drives the 
computed phase to zero. 

One aim of this study was to look for center-to-limb variations in several oscillation 
parameters, and then to express these as a function of height through the sola r atmosphere. 



Some studies of oscillations have been d one at different limb positions before (jSolanki et al. 



19961 : lUitenbroek. Noyes fc Rabinlll994j ) however having the entire disk available in this data 
set provided a considerable advantage. Since the data contain about 60,000 pixels across the 
solar disk, binning in annuli of A/i = 0.01 enabled us to average our oscillation parameters 
using between 600 and 1200 pixels. This substantially increases the signal-to-noise ratio in 
our data when compared with a single measurement at a given value of fi. 

Finally the two-dimensional aspect of this data enabled the production of diagnostic 
diagrams (or 1-z/ diagrams). Many studies have been done to examine the I-V phase us- 



ing d i agnostic diagrams with other spectral lin es (lAl. Bendlin fc Kneen Il998l : lOliviero et al. 



I999I : IStrauss et al.l Il999l : iDeubner et al.l Il992l ) but our data set sampled a unique atmo- 
spheric height over a unique range of spatial wavelengths. For this analysis, a section of 
roughly 1000x1000 arcseconds from the disk center was extracted from both the velocity 
and line depth time series. No correction was made for projection effects or differential solar 
rotation, nevertheless valuable information regarding the spatial and spectral dispersion of 
the observed oscillations is obtained from this analysis. 



4. Discussion 

The shifts of the CO line centers across the solar disk are consistent with the expected 
solar photospheric rotation velocity. After removing solar rotation, a low pass temporal 
filter was applied to the time series in order to extract the supergranulation component of 
the disk velocity. In these images the supergranulation cells are clearly visible, and a fit 
to the outfiow velocities (both positive and negative lobes) suggests a horizontal outflow 
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velocity of 186 m s~^. This is in agreement with the distribution of outflow velocities seen 
from Hinode feature tracking measurements in the UV continuum channels (Tian et al 2010) 
which originate from the same approximate heights. The behavior of the supergranulation 
outflow Doppler signals with fi suggests a height variation in the outflow speed, which will 
be investigated in future work. An analysis of the line bisectors of the CO lines showed no 
variations, which was likely due to the rather low spectral resolution achieved with these 
observations. 



4.1. Spatial Variations 



Using the low pass temporal filter, maps of the line depth across the solar disk were made. 
The data show that on a global scale the CO line absorption appears weakened in magnetic 
features associated with the solar magnetic network; confirming results from previous higher 
resolution studies ( lUitenbroek. Noyes fc Rabinlll994j : lAyres fc Rabinlll996l ). In Figure 2 we 
show an image produced frorn a ma gnetogram taken on the same day from the NSO SOLIS 
instrument (i.e. iJones et al.l (|2002[ )): here the absolute value of the magnetogram is taken. 



and then the map is smeared with box-car averaging to a resolution of about 10 arcseconds (a 
little sharper than the expected CO map resolution). In Figure 2 we plot the mean CO line 
depth (averaged over the three spectral lines) across the solar disk, excluding pixels near the 
limb where seeing and fitting noise dominate the map. There is a very clear correspondence 
between the position of stronger magnetic fields in the SOLIS map with the regions of weaker 
CO line depth. 

The beating of solar p-modes with slightly different frequencies complicate the simple 
cross-correlation of the line depth and velocity time series. This is especially true if the 
power distribution among the oscillation modes varies in the line depth compared to the 
velocity measurements, and such a different power distribution is a natural consequence 
of non-adiabatic atmospheric conditions. Furthermore, if the I-V phase varies with oscil- 
lation frequency we expect even more difficulty with a simple cross-correlation technique. 
Nevertheless computing the cross-correlation was very robust. As previously noted, the 
cross-correlation lag was measured in units of seconds of time, but then was converted to an 
angular phase shift assuming a dominant oscillation period of five minutes. In most cases, 
the peak of the cross-correlation function occurred with the line depth leading the velocity 
by about 90 degrees; in some cases however the peak in the cross-correlation function at -270 
degrees was higher. To unwrap the measurements, we forced our analysis software to fit the 
peak near 90 degrees. Figure 2 shows a map of the phase shift across much of the solar disk. 
As the Sun had very low activity, only small variations are noted; however these variations 
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are significant, and are not tightly correlated with the line depth map of Figure 2b. The 
phase map shows something different than the line depth map, and presumably the varia- 
tions represent regions of more or less adiabatic response in the solar atmosphere. As we 
will discuss in detail later, one interpretation of this is that the phase map show regions with 
faster and slower thermal relaxation. 

The Fourier analysis of the line depth and velocity time series was done to more fully 
investigate the properties of the oscillations. Both line depth and velocity showed power 
spectrum peaks at the roughly 3.3 mHz, equivalent to a period of five minutes. Unlike other 
studies of oscillations with these spectral lines, these data show no power for oscillations 
near 5 mHz (3 minute period). It must be stated that the signal-to- noise ratio of these data 
benefit significantly from the ability to average across many pixels from the solar surface; 
in some of the following discussions we average over tens of thousands of pixels, and so our 
power spectra has a factor of 100 times better signal-to-noise than previous works which 
report spectra from single positions on the solar disk. 

The oscillation power integrated from 3.0 < u < 4.2 mHz using the CO velocity power 
spectrum measured pixel-by-pixel across the solar disk is also shown in Figure 2. Regions of 
high power and lower power are clearly seen. By examining the 500 weakest and strongest 
line depth pixels on the CO line depth map, (also corresponding to regions of stronger and 
weaker line-of-sight magnetic fields) a power deficit in magnetic regions is seen as the p-mode 
power (after background subtraction). In the regions of weaker line depth pixels (stronger 
magnetic fields) the power is about a factor of 0.5 times the power in the strongest line 
depth pixels. There is also a weak correspondence between oscillation power and phase, 
with regions of high (low) power showing non-adiabatic (adiabatic) phases. 

The interest in observing the full solar disk with these observations was to look for both 
large and small spatial scale variations of the oscillations parameters. The variation from 
center-to-limb can be used to characterize the height variation within the solar atmosphere 
by using the Eddington-Barbier relationship. Here the approximation is made that the 
source function arises solely from a thin layer at r = 1, and that the height of this layer 
changes across the solar surface with fi. Using our measured I-V phases, it becomes possible 
to convert from r to height z in the solar atmosphere. As discussed by Kopp et al. (1992), 
the formation of the infrared continuum in the solar atmosphere is more easily understood 
than the formation of spectral lines, as the continuum is formed in LTE and the dominant 
source of opac ity is H~ free-free absorpti on. Given a solar atmospheric model, (they use 



model C from IVernazza. Avrett fc Loeserl (jl98ll )) they determine the formation heights of 



their 50, 100, 200, and 400 fi IR continuum channels as 340, 420, 480, and 600 km above 
z = 0. By subtracting their observed phase shifts from 180 degrees (since they use Doppler 
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velocity rather than atmospheric velocity), we can compute I-V phase lags for these heights 
as 122, 122, 109 and 95 degrees respectively (with errors of about ±10 degrees). 

Figure 3 shows the mean I-V phase shift as a function of fi from our CO observations. 
Here we have averaged the phase spectra from 3.0 to 4.0 mHz, and averaged over annuli on 
the solar disk with a width of A/i = 0.01. The I-V phase drops from about 116 degrees at 
disk center to about 100 degrees at /i = 0.5. These values are clearly within the range of 
phases observed by Kopp et al. (1992) and we can now determine the atmospheric heights 
from which they arise, assuming that equal phases in our data correspond to equal heights 
in the IR continuum data. In detail, we use the Eddington-Barbier relation to replace /n(r) 



with I njp,). We determine the relationship between z and (pi-v from the data from iKopp 



( 119901 ). and then use the relationship between and (j)i-v in our observations to determine 
a ln{fi) vs z function for the formation height of our CO line depth. Our CO I-V phases 
between 0.5 < /i < 1.0 correspond to heights between 560 > z > 425 km. and the scale 
height which we derive for the variation of the CO I-V phase gives Hi = 195 km. 

We can now examine the power in the CO oscillations as a function of disk position, and 
therefore height. Figure 4 shows the CO velocity power spectra averaged over the spatial 
positions from 0.99 < fi < 1.00. It shows a significant background power which varies 
with frequency, and a broad peak near 3.3 mHz which represents the well-known p-mode 
oscillations. We fit a polynomial function to the log of this background using the frequency 
ranges of 1.5 < u < 2.0 mHz and 5.0 < u < 8.0 mHz and interpolate the background 
power through the other frequencies. After subtracting this background, we fit a Gaussian 
function to the power spectrum. The fits to the power spectrum are shown as the solid 
line in Figure 4. The fits to the oscillation power across the solar disk were examined. The 
peak of the p-mode power given by the Gaussian fit was corrected for two projection effects: 
first, the power was divided by fi to compute the oscillation power assuming purely radial 
motions, and second the power was corrected for the smearing in spatial wavenumber caused 
by projection effects as one moves from the center to the limb. (This correction, basically 
a distribution of power with spatial wavenumber, was computed by averaging the velocity 
signal in pixels near disk center and then computing the power spectra of these averaged 
velocities.) After making these two correction, we transform the fi coordinate into a height 
coordinate z and examine the power changes with height. 

Figure 4 shows that the total power (p-mode and background) increases with height, 
whereas the p-mode power al one decreases with height. The increase of the total power is con- 



sisten t with previous studies ( iRuiz Cobo. Rodriguez Hidalgo fc Colladoslll997l : ISimoniello et al 



20081 ) and supports the idea that the velocity oscillation power must increase with height 



as the density drops. An exponential fit to this background power gives a scale height of 



- 11 - 



H2 = 230 km, consistent with values used for chromospheric models. By examining the 
p-mode power alone, after subtracting the background power, we find that it decreases with 
height, as would be expected from an evanescent solution to the wave equation. This de- 
crease is very rapid; the scale height is much smaller and consistent with photospheric scale 
heights, = 90 km. 

Another parameter derived from the power spectra fitting is the peak central frequency 
for the p-mode oscillations. While many studies find 3-minute period oscillations at chro- 
mospheric heights and lower (using these CO spectral lines) we find no significant power 
peaks at that frequency. Moreover, the central frequency position of the fit to the p-mode 
power does not change with height through the range of 560 > z > 425 km. No 3-minute 
period power is seen in an analysis of the line depth oscillation power spectra, and the central 
frequency of the p-mode oscillations as measured with the line depth also remains constant 
with disk position. 



4.2. Diagnostic Diagrams 

In Figure 5 we present diagnostic diagrams, or 1-z/ diagrams of several parameters from 
the CO data; these diagrams are produced using the three dimensional Fourier transform of 
the region near the center of the solar disk as discussed in Section 3. A simple correction for 
solar rotation was made by aligning the ridge structure for positive and negative frequencies 
before coadding; but no correction is made for projection effects. The CO velocity l-u diagram 
clearly shows several ridges of power, and these six ridges are at the position of the well- 
known n=2 through n=7 ridges for the global oscillation modes. The CO line depth diagram 
is noisier, but shows a couple of ridges at the positions of the global oscillation modes. In 
both parameters, no distinct ridge structure is seen at high frequencies. At 4 < z/ < 5mHz 
the ridges fade into the background power; and no other structure is seen at the u = 5mHz 
frequency where some studies have reported to have seen chromospheric oscillation power. 

The coherence between the velocity and line depth in Figure 5 also shows strong ridge 
structure. While the resolutions Al and Au not especially high, there is some evidence that 
the coherence betw e en th e global p-mode ridges does not go to zero. This is consistent 



with lOliviero et al.l (119991 ) who show significant background power at the inter-ridge posi- 
tions in measurements of global p-modes near these spatial frequencies. The coherence also 
shows significant non-zero values at high temporal frequencies, up to about u = 6.5mHz, 
especially at low spatial degree. Although there are no ridges or peaks at these frequencies 
suggesting distinct oscillation modes, this large coherence value suggests that waves at this 
chromospheric oscillation frequency may be seen with these CO lines. 
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The I-V phase diagnostic diagram in Figure 5 shows a correlation with the coherence 
diagram: the phase values are noisy where the coherence is small, and the phase values change 
smoothly where the coherence is large. There is very little evidence for a ridge structure in 
this data, and the phase is rather constant at all spatial wavelengths. The phase does show a 
gradual decrease as one moves to lower temporal frequencies, from about = 130 degrees at 
u = 2.5 mHz to roughly (j) = 90 degrees at i' = 4.5 mHz. (It should be noted that trend this 
agrees with 180 — from lOliviero et al.l (119991 ). assuming that the phase in that work was 
measured using Doppler velocities.) There is no evidence for a sh arp change in the ph ase 
or a plateau of low phase below u = 2 mHz as seen in other work ( iDeubner et al.lll990l ): as 
suggested by the low coherence at these frequencies, this may simply reflect that the phase 
in these observations is not well measured here. 



4.3. Phase Spectra Fits with Analytic Model 



The Kopp model (IKoppI Il990l : iKopp et al.l Il992l ) is based upon several simplifying as- 
sumptions. First it uses hydrostatics, thus ignoring the gas dynamics and the magnetic 
fields which are present in the solar temperature minimum and chromosphere (effectively 
the magnetic field is assumed to be less than about 200 Gauss). The model is based upon 
an isothermal atmosphere which simplifies the temperature profile in this region of the at- 
mosphere. And finally, the model assumes a plane-parallel atmosphere which may have 
implications with these full-disk me asurements. Ne vertheless, the Kopp model is one of the 
few analytic models (in addition to IWorralll ( 12002| )) to date which describes the I-V phase 
for solar oscillations in physical parameters; because of this it is still very useful. 



KoppI (jl990l ) develop expressions for the I-V phase in terms of the atmospheric acoustic 



cutoff frequency Uc, the scale height H the ratio of specific heats 7 and the thermal relaxation 
frequency Ur- The work makes no assumption regarding the actual physical mechanism for 
the thermal relaxation process, it just assumes that the p r ocess follows a Newton cooling 
law. What is explicitly lacking from the analysis of iKoppI (Il990[ ) is a direct expression for 
the I-V pha se spectrum, 0(i^). In the Appendix below, we continue the development from 
KoppI ( I1990I ) one extra step to derive and expression for 0(z/) which depends solely on the 
physical variables listed above. Our expression is listed in Equation 5 of the Appendix. 

We use our expression for 0(z/) to fit our I-V phase spectra as determined at different disk 
positions from our data. The points Figure 6 shows a sample I-V spectrum measured near 
disk center. An adiabatic model with no thermal relaxation would produce a flat spectrum 
with a 90 degree phase shift at all frequencies; this is clearly contradicted by the data which 
show a decrease from about 120 degrees to about 105 degrees from frequencies of about 2.8 to 
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4.0 mHz (these frequencies are where the I-V coherence is high). Following iKoppI ( 119901 ) we 
can assume that the acoustic cutoff frequency is C(;c/(27r) = 4.5mHz and do a one parameter 
fit to the data to determine the best value of the thermal relaxation frequency Ur- This Id 
fit is shown as the dashed line in Figure 6; it has a high value and does not produce a 
satisfying fit to the data points. If we let the acoustic cutoff frequency be a free parameter 
and do a 2d fit to the spectrum, we significantly reduce the value of and derive a more 
satisfying fit to the data points. 

An examination of the surface which results from a 2d fit using c^c and Ur reveals 
that in some cases the analytic function only weakly constrains Uc- This combined with 
the inherent non-linearity of the atari function have foiled our attempts to fit this data 
using standard non-linear least squares software. Instead, for each I-V spectra at a given 
disk position (fi) we simply examine the over the range of < Ur/{2TT) < lOmHz and 
< Uc/ (27r) < lOmHz and find the point of minimum x^. We estimate the 68% confidence 
ranges by u sing specific i ntervals on this surface according to standard numerical analysis 
techniques (jPress et al.lll994j ). 

Figure 6 shows the values we derive from fitting the I-V phase spectra over the range of 
0.25 < fi < l.OmHz. As in the case of Figure 6, we plot the results for the Id fit using the 
dashed line, and the 2d fit results using solid lines. In the Id case we see that the thermal 
relaxation frequency Ur drops from the center of the disk to the limb; this suggests that as 
the data probe higher layers of the solar atmosphere the thermal relaxation process becomes 
less important and the atmosphere becomes more adiabatic. The same behavior for Ur is 
shown with the 2d fits where the relaxation frequency decreases at /x decreases, although 
the values are consistently O.SmHz higher than in the Id fits. And finally, the value for 
the acoustic cutoff frequency is around 4.5mHz, although the data suggest the value may 
significantly drop near the solar limb (and thus at higher heights in the solar atmosphere). 

The values of 0.5 < Ur/{2TT) < 2.8mHz correspond to a cooling time scales between 
2000 and 360 seconds which are longer than the values of 290 to 76 seconds fo r this range 



of hei ghts in the solar atmosphere measured with differe nt fitting methods in iKopp et al. 
( 119921 ). While the relaxation rate (defined by iKoppI ( 119901 )) of ^ which we obtain is similar 
to the values determined by iKoppI (119901 ) , our values of cOr are somewhat lower since our 
fitted values of Ur are somewhat lower than the 4.5mHz assumed in their work. 
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5. Appendix: Analytic Expression for the Phase Spectrum Including Thermal 

Relctxation 

Kopp (1990) develops an analytic expression for the dispersion relation for a thermally 
relaxing chromosphere. Assumptions in his simple model include an isothermal atmosphere 
near the temperature minimum, no magnetic fields, the plane-parallel approximation, and a 
Newton's cooling law for the thermal relaxation process. While these assumptions may be 
a poor reflection of the actual physics of the solar chromosphere the work gives one of the 
few clear analytic expressions for the I-V phase relation published to date. We continue the 
discussion from Chapter 5 of Kopp (1990) a few steps to determine an expression for the I-V 
phase as a function of frequency from his model. 

Kopp starts with the conservation of mass, the equation of motion, a thermal equilib- 
rium equation including Newton cooling, and the ideal gas law (see his eq. 5.4 through 
5.7). He linearizes the equations aroimd an equilibrium solution, assumes vertical per- 
turbations in pressure, density and velocity, (where the vertical velocity perturbations are 
V = wz'exp**^'^^"'^^^-* with w as a unitlcss constant) and determines the dispersion relation for 
vertical wavenumber kz as a function of oscillation frequency cu as follows (Kopp's equation 
5.25): 

ife.^^f-l±Jl.- ,^-.f-/"\ (^X] (1) 

Cs \ V 1 ~ l{uJr/l(^) \^cj J 

where Uc = ^ is the acoustic cutoff frequency, Cg is the sound speed, 7 is the ratio of heat 
capacities, H is the scale height and cu^ is the thermal relaxation frequency. 

By assuming that can be written in the form 

kz — ky -\- iki (2) 

where kr and ki are purely real values, and by writing an expression for temperature per- 
turbations in a similar form to the velocity perturbations 5T — ^Tq Q-^jp^i'^^-^zz) ^j-^-j^ ^ 
a unitless constant, Kopp relates the temperature and velocity perturbations as (equation 
5-29): 

(3) 



-fujkr/cl _ ■ ( (l^l^^s)iki - l/H) ki 



Finally, the I-V phase can be determined by (equation 5.30): 



A;2 + [ki - l/Hf u J \kl + {k, - l/Hf + uj 



= arctan ( ) (4) 



To develop an explicit expression for the phase as a function of frequency, we define 
the ratios a = ui/uic and h = our/ujc and gather terms in Equation 1. Next we reduce the 
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square- root of the complex number and solve for kr and ki, keeping only real solutions for 
these two coefficients. Finally we substitute the wavenumbers into the phase expression to 
derive: 

ra\h - II H) - h{kl + {h - l/Hf) 



(f) — arctan 



-7 , 



(2H)2 krd^ kr(^k'^ -\- {ki \j Hy^^ 



(5) 



where: 



ki 

kf 



1 
2H 
1 



(1-/) 

1 X2, 



f 



X2 + 



•^2 "I" -^3 



2 2 7 ' 

Xi = a 7 + 

,2„,2/i „2\ I 1,2/1 „2 



and 



X2 = aV(l-«') + &'(l-«'7) 



X2 = a^&7(7 - 1) 



As a check for this solution, we reproduce Figure 5.2 of Kopp (1990) using values for 7 = 5/3 
and H — 150 km. We also select the solution for 'upwardly locahzed waves' (not shown 
here) and successfully reproduce Kopp's Figure 5.3. Section 4.3 of this paper discusses the 
agreement between our Equation 5 and our observations and the detailed fitting we do using 
this expression to derive the behavior of cUc and cUr in the solar atmosphere. 
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Fig. 1. — A spectral frame from the solar disk scan at 1647UT on 16 Sep 2009 using the 
NAC. The three CO lines used in this study are labeled; the other features in the spectrum 
are all C O lines except for the broad atmospheric absorption at 4666. 8nm. (Identifications 
are from [Wallace fc Livingston! ( l2003l )). The spectrograph slit roughly lies along the central 
meridian of the solar disk in this frame. 
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Fig. 2. — Maps of the absolute values of the SOLIS magnetogram (smeared to lower spatial 
resolution), the CO line depth, the CO Doppler velocity p-mode power, and the CO I-V 
phase compute from the cross-correlation of the line depth and atmospheric velocity time 
series. Solar north is up and solar west is to the right in each of the images. Each image 
has circles at 0.2, 0.4, 0.6 and 0.8 solar radii, and the center of the disk is near the left 
edge of each image. (The power and phase maps are only displayed up to 0.8 solar radii.) 
The two smaller circles show regions of reduced CO line depth, which correspond to regions 
of vertical magnetic field in the SOLIS magnetogram. Both regions show reduced p-mode 
power, but while the phase lag in these regions is weakly correlated (showing phases close to 
90 degrees) it also shows uncorrelated structure. The phase map is scaled from about 130 
degrees (white) to 90 degrees (black). 
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Fig. 3. — The I-V phase shift averaged from 3.0 < v < 3.4 mHz as a function of position 
across the solar disk, from the center of the disk /j, — 1 toward the hmb. The value of (j)i-v 
would be equal to 90 degrees if the solar atmosphere responded adiabatically, but this is 
clearly not the case over much of the solar surface. Using the Eddington-Barbier relationship, 
we can connect a decreasing value of /i with larger heights in the solar atmosphere; and so 
the trend of to move towards 90 degrees as /i decreases means that the atmosphere 
behaves more adiabatically with increasing height. 



- 22 - 




Fig. 4. — The disk center (// > 0.99) CO Doppler shift shows a strong peak in the power 
spectrum (left) at u — 3.2 mHz, consistent with the solar p-mode frequencies; the symbols 
show the measured power. The background power is fit with a simple polynomial and drops 
smoothly with increasing frequency, and the p-mode peak can be roughly fit with a Gaussian 
function. No evidence is seen for 3-minute period oscillation power near z/ = 5 mHz. The 
figure on the right shows how the total power (background plus p-mode peak) and the power 
in the p-mode peak behave at different heights in the solar atmosphere. The power spectra fit 
parameters from different disk positions were corrected for projection and smearing effects, 
and then the values were used to determine a physical height corresponding to a 

particular value for /i (see text for details). The total power increases with height with a 
scale height oi H — 230 km, and the p-mode peak power drops with height, consistent with 
an evanescent wave and with a scale height of only if = 90 km. 
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Fig. 5. — Diagnostic l-u diagrams for various parameters in this study, taken from a region 
near the center of the solar disk. From left to right we present the CO velocity diagram, 
the CO line depth diagram, the coherence between the line depth and velocity, and the 
phase between the hne depth and velocity. The global p-mode ridges from measurements of 
the NSO GONG project for n—2 through n—7 are displayed velocity diagram, and exactly 
correspond to power ridges seen there. Neither the velocity nor the line depth diagrams show 
significant power near i/ — 5 mHz, although the coherence diagram has a significant value 
above the noise up to about i/ = 6 mHz, and the phase shows signal levels above the noise up 
to these higher frequencies as well. The values for the phase start near = 130 degrees 

at low frequencies u = 2.5 mHz, and drop to values of about 4>i-v = 90 degrees at the higher 
frequencies, = 4.5 mHz, 
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Fig. 6. — Fits to the phase spectra are shown in the figure on the left, and the changes 
of the fit parameters across the solar disk are plotted in figure on the right. The phase 
spectra 01-1/ (z/) were fit with an expression developed in the Appendix in two ways: a "1- 
dimensional" fit was done with a constant acoustic cutoff frequency, ujcj (27r) = 4.5 mHz and 
finding the minimum ^ value to determine the thermal relaxation frequency, w,.. This fit 
is shown as the dashed line, and has a large value of . The behavior of the fit parameter 
LOrj (27r) across the solar disk from these " Id" fits is shown as the dashed line in the figure 
on the right. A " 2-dimensional" fit was done using a surface to fit both the acoustic and 
relaxation frequencies. This provides a better fit to the phase spectrum on the left, and the 
fit parameters are shown as they vary across the solar disk (and presumably with height in 
the solar atmosphere) on the right. In the right figure the error bars are 68% confidence 
levels determined from the variations of . 



